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Z2 spin liquid and chiral antiferromagnetic phase in Hubbard model on the honeycomb 
lattice: duality between Schwinger-fermion and Schwinger-boson representations 

Yuan-Ming Lu 1 and Ying Ran 1 

department of Physics, Boston College, Chestnut Hill, MA 02467, USA 
(Dated: July 19, 2011) 

In our previous work-, we identify the Sublattice-Pairing State (SPS) in Schwinger-fermion rep- 
resentation as the spin liquid phase discovered in recent numerical study on a honeycomb lattice . 
In this paper, we show that SPS is identical to the zero-flux Z2 spin liquid in Schwinger-boson rep- 
resentation found by Wang^ by an explicit duality transformation. SPS is connected to an unusual 
antiferromagnetic ordered phase, which we term as chiral-antiferromagnetic (CAF) phase, by an 
0(4) critical point. CAF phase breaks the SU(2) spin rotation symmetry completely and has three 
Goldstone modes. Our results indicate that there is likely a hidden phase transition between CAF 
phase and simple AF phase at large U/t. We propose numerical measurements to reveal the CAF 
phase and the hidden phase transition. 



PACS numbers: 

Introduction In a recent numerical study 2 , it is 
found that the Hubbard model on the honeycomb lattice 
hosts a spin disordered insulating phase in the neigh- 
borhood of the Mott transition, which does not break 
any physical symmetry. If this is the case, this phase 
should be a novel spin liquid phase with fractionalized 
excitations 4 . Similar fractionalized spin liquids have been 
shown to exist in various artificial models^—, but so far 
there is no simple and hopefully realizable Hamiltonian 
that hosts such exotic phases. This remarkable numeri- 
cal study makes an important footstep along finding spin 
liquid phases in correlated electron systems. 

There are a lot of different spin liquids on the hon- 
eycomb lattice, characterized by different topological or- 
ders, or different Projective Symmetry Groups (PSG) 11 . 
For example, in a previous study we show that there can 
be 128 distinct spin liquid phases within the Schwinger- 
fermion representation 1 . Which one is realized in the 
Hubbard model? In the numerical study^, it is shown 
that the spin liquid phase have a full energy gap, and 
is smoothly connected {i.e. through a continuous phase 
transition) to both the semi-metal phase for small U/t 
and the Neel phase for large U/t. These three conditions 
strongly restricts the candidate spin liquid phases. 

In our previous study by Schwinger-fermion approach^, 
we use only two of the three conditions and show that 
there is only one natural spin liquid, coined the Sublattice 
Pairing State(SPS), which has a full energy gap and can 
be smoothly connected to the semi-metal phase. Is SPS 
compatible with the third condition? In other words, can 
SPS be connected to a magnetic ordered phase by a con- 
tinuous phase transition? Describing ordinary magnetic 
ordered phase in Schwinger-fermion approach has been a 
puzzle for a long time. In this work we show a solution, 
which allows us to study phase transitions between SPS 
and a magnetic ordered phase. 

SPS is a fully-gapped Z-i spin liquid on the honeycomb 
lattice. Its mean-field fermionic spinon band structure, 
after choosing a proper gauge, is given as follows: (see 



FIG. W) 

H SPS = *1 2^ fiafja +h 2^ fiafjoL 

-ME /-•/» + A E ^Itfl+h.c. (1) 

where t\,% are real numbers. In Schwinger-fermion 
approach, /-spinons are coupled to an SU{2) gauge 
fiel d 11 ' 12 . However due to non-zero £2 and A, this 517(2) 
gauge degree of freedom is reduced to Zi through Higgs 
mechanism. Thus at low energy /-spinons in Eq.([T]) arc 
coupled to a dynamical Z2 gauge field and stay in the 
deconfincd phase. 

Recently Wang 3 identified the 0-flux state as the most 
promising spin liquid phase in the Schwinger-boson ap- 
proach. 0-flux state is also a Z^ state with a full energy 
gap, smoothly connected to an antiferromagnetic phase 
through an 0(4) critical point. However it is not clear 
whether the 0-flux state can be smoothly connected to 
the semi-metal phase. Can be SPS related to the 0-flux 
state? The two states are described in two completely dif- 
ferent languages: one in Schwinger-fermion and the other 
in Schwinger-boson, and the relation between these two 
representations has been a longstanding problem. In this 
paper, we find that strikingly, SPS and 0-flux state are 
identical by an explicit duality transformation in the low 
energy effective theory. 

We also find out that the antiferromagnetic phase con- 
nected to the 0-flux state (or SPS) is rather unusual 
and not the simple Neel phase, because it breaks the 
SU(2) spin-rotation symmetry completely and has three 
Goldstone modes. On the other hand the magnetic 
order is still colinear. We dub this phase as chiral- 
antifcrromagnetic(CAF) phase. In CAF phase, aside 
from the usual antiferromagnetic spin order parameter 
N = (—l) ls Si where i s — 0,1 for A and B sublattices 
respectively, there is another vector-chirality spin order 
parameter n = 5Z<<«>> VijSi x Sj whose expectation 



value satisfies (n) _L (N), and v± 
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FIG. 1: Proposed schematic phase diagram of the Hubbard 
model on the honeycomb lattice. The hidden phase transition 
might happen at a large U/t, which is not reached in the 
numerical study^. 



a left (right) turn when going from site j to i as shown by 
the arrows in FIG. |3Jd. Since the usual AF phase should 
exist in the large U/t limits, our results suggest a hidden 
phase transition, which might happen in the "Neel" or- 
dered phase of the numerical study 2 - or at larger U/t not 
studied before. In other words the numerical study may 
not distinguish CAF and simple AF phases. We propose 
the schematic phase diagram as shown in Fig[TJ 

Continuous phase transition from SPS to CAF 
phase We first discuss the continuous phase transition 
from SPS to CAF phase. Describing an ordinary mag- 
netic ordered phase in Schwinger-fermion approach is 
highly non-trivial: this is because a Schwinger-fermion 
mean-field ansatz has at least an unbroken Z2 gauge 
symmetry and to describe a regular magnetic phase, the 
gauge degree of freedoms must be confined. The first 
demonstration of a regular magnetic-ordered phase in 
Schwinger-fermion approach is given in Ref. 14. It is 
shown that the easy-plane antiferromagnetic order (XY 
order) on the honeycomb lattice is described by a quan- 
tum spin Hall (QSH) band structure of the fermionic 
spinon along S z direction coupled with a dynamical U(l) 
gauge field, where f {S z = ^) spinon band has a Chern- 
number v = 1 while \. (S z = — ^) spinon band has a 
Chern- number v = — 1. Because of the QSH effect, the 
gauge fluctuation is bound to S z spin density fluctua- 
tion, and the Goldstone mode of the easy-plane Neel or- 
der is nothing but the photon of the U(l) gauge field. 
The long-range spin-spin correlation function in the Neel 
phase is dual to the long-range monopole-monopole cor- 
relation function in the Coulomb phase of the U(l) gauge 
field. Indeed one can compute the monopole quantum 
number— to show the spin order pattern is antiferromag- 
netic. 

Armed with this result, we now consider an SU(2) spin 
rotation symmetric system. CAF phase is described by 
spinon band struture in the presence of a fluctuating 0(3) 
QSH order parameter n, coupled with a U(l) gauge field. 
Its mean-field ansatz is: (see FIG. [3})) 



H CAF — *1 7 y fiafjet +t2 2^ fiafjc 

-u^fLfit+n- J2 WijfLZotpfjp, (2) 



«IJ» 



Note that the non-vanishing ti term, as well as the QSH 
order ft, breaks the SU(2) gauge symmetry down to f (1). 




FIG. 2: The honeycomb lattice and its Bravais lattice vector 
5i,2- rfi,2,3 are the three vectors defined in the effective La- 
grangian of z boson for the SPS-CAF phase transition. Two 
generators of symmetry group are also shown: 60° rotation 
C% around the center of a plaquette and the horizontal mirror 
reflection a. 



In the Dirac limit (£2, \n\ = tn -C t\) we can write down 
the effective Lagrangian of CAF phase in imaginary time: 

C =^70(0^ - m M )7 M V + mn- ip^Mip 



>f, 



2-iy-v 



2u 



(dph) 



(3) 



where ^ is a 8-component complex fermion describing the 
long wavelength part of /-spinon around the two Dirac 
points K and K' = —K in the Brillouin Zone: ipk = 
(fK+k,A-, fK+k,B, fx'+k.B, fK'+k.A) and we neglect spin 
indices. 7 M are Pauli matrices in the sublattice space: 

7o = M3>7i = M2,72 = -Mi- M = 70^ are the QSH 
masses, and space/time are rescaled so that the Fermi 
velocity is one. Non-universal coupling u describes the 
fluctuation of the QSH order parameter: unit vector n. 

There are three gapless modes in CAF phase: two h 
fluctuating modes and one photon mode. The photon 
mode is in-plane spin wave of anti-ferromagnetic order 
N (N _L n), and the spin SU{2) symmetry is completely 
broken. The combination C§ ■ T of the 60° space ro- 
tation around the hexagon center and time-reversal T 
leaves both order parameters invariant. This symmetry 
indicates that the magnetic order in CAF phase is still 
collinear. 

Comparing Eq.® with Eq.([T|), s-wave pairing A of 
spinons in SPS phase is replaced by the 0(3) QSH order n 
in CAF phase. If we group these orders together into a 5- 
component vector V = (ReA,ImA,n) and ignore gauge 
fields for the moment, as pointed out by Grover and 
Scnthil 15 , fluctuations of V has a Wess-Zumino-Wittcn 
(WZW) term after integrating out fermions. This WZW 
term allows us to construct an 0(4) phase transition be- 
tween SPS and CAF. 

WZW term is a topological Berry phase in the non- 
linear sigma model of V = («i, «2, V3, V4, V5) (\V\ = 1) in 
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FIG. 3: (color online) Mean-field ansatz of (a) SPS phase and 
(b) CAF phase in terms of /-fermion. ti,2 are real fermion- 
hopping parameters, while s-wave pairing order parameter A 
can be complex. Vector n is the QSH order parameter in the 
CAF phase. ny — 1 if * — ¥ j is along the arrow direction. 



2+1 dimension 16 : 



5' 



2?ri 



wzw 



Area{S A ) J 



dp j d 3 xe abc deV a dpV b dtv c d x VddyV e 

(4) 



where Area(S 4 ) — ^-, and V(x 7 y,t 7 p) is a smooth ex- 
tension of space-time configuration V(x,y,t) to 4-disk: 

V(%, U, t, P = 1) = V(x, y, t) and V(x, y, t, p = 0) = V is 
a fixed vector. The physical meaning of the WZW term 
is that a Skyrmion (anti-Skyrmion) of h in two space 
dimensions actually carries fermion charge 2 (—2). And 
the hedgehog instanton of n in 2+1 dimension creates a 
charge-2 s-wave fermion pair. Due to this WZW term, 
as shown in Ref. Ila . a direct phase transition between 
a QSH insulator and an s-wave superconductor on the 
honeycomb lattice becomes possible. 

Let us keep the WZW term in mind and study the 
CAF to SPS phase transition. Starting from the effec- 
tive theory of the CAF phase Eq.Q, it is convenient to 
introduce the CP 1 representation of the h order param- 
eter: n = w'aw where w — (wi,w 2 ) T are two complex 



numbers satisfying |tui| 



w 2 



1. This representation 



has a 17(1) gauge redundancy and thus w-bosons couple 
to a U(l) gauge field A^ at low energy. After integrating 
out the /-spinon of Eq.© the effective Lagrangian is: 



L =|(0„ - iA,)w\ 2 + rM 2 + S \w\ A + -L/*, 

Z 9a 

+ tt^tF -\ — euv\Aud u a\ 

2 9 A n 



(5) 



where the constraint |u>i| + |u>2| = 1 1S softened . 
ffj,u = dfj,a„ - <9„a M and F^ v = d^A v - d v A^ are the field 
strengths of the two gauge fields. A key observation of 
the present work is the mutual Chern-Simons (CS) term 
between A^ and a^. While an explicit derivation of this 
term is given in Appendix [Al it can be easily understood: 



a Skyrmion of n is well-known to be represented as a 27r 
Ap gauge flux. Now the WZW term dictates this 2tt A^ 
flux to carry charge-2 of a M gauge field, which is exactly 
described by the mutual CS term. 

What are the phases described by the effective La- 
grangian Eq.([5])? When r < 0, w-boson condenses and n 
is ordered, so the system is in the CAF phase. Note that 
in this phase the mutual Chern-Simons term does not 
qualitatively modify the low energy gauge dynamics: the 
A^ gauge fields are gapped out due to the Higgs mech- 
anism. Remarkably, when r > 0, w-bosons are gapped 
and h is disordered, and the system is actually in the Z 2 
SPS spin liquid phase. The U(l) mutual Chern-Simons 
term or the WZW term is essential to make this iden- 
tification. The relation between mutual CS theory and 
Z 2 gauge theory was firstly discussed in Ref. [17|. Here, 
based on the WZW term and monopole physics, we are 
able to further identify the PSG of this Z 2 state. 

In the n disordered phase, the mutual CS term opens 
up mass gaps for both A^ and a M gauge fluctuations^: 
TO a ~ 9A> m A ~ 9a- Meanwhile, the real time 
(Minkowski spacetime) equations of motion (with source 
term —a^j^ — A^J^) become 

U = ^d"U + -e a ^d a A^ 

9a * 

Jp = \d v F vll + -e a ^d a a + iw^w, (6) 
g A * n 

where the last term is the w-boson current i[w'(dij, — 
iA^)w — w(d ll + iAn)w']. These equations of motion in- 
dicate that one flux quanta of a^(A^) carries two units 
of A^icif,,) gauge charge. And one unit of A^ 1 (a fl ) gauge 
charge in turn carries 7r-flux of a tl (A^ t ) gauge field. As 
a result, /-spinon sees the w-boson as a 7r-flux and vice 
versa. These are the two gapped fundamental excitations 
in a Z 2 gauge theory: gauge charge and vison, and they 
are dual degrees of freedom. If we define the /-spinon (w- 
boson) as the Z 2 gauge charge, w-boson(/-spinon) will 
be the vison. From the viewpoint of /-spinons, we just 
showed that their visons, i.e. w-bosons carry spin quan- 
tum number. 

Now let's discuss the monopolcs of a M and A^ in Eq.©. 
In n ordered phase (CAF phase), where A^ is Higgsed 
out, the a M monopole is bound with a spin flip due to the 
QSH effect. Therefore a^ monopole events are suppressed 
in the CAF phase by spin rotation symmetry and a M stays 
in Coulomb phase. However in the fi disordered phase, 
there is no QSH effect and monopole events of both a M 
and A^ are allowed. Let us denote their monopole cre- 
ation operators as V^ and V\ respectively. What are the 
physical consequences of these monopole events? 

By equations of motion it is shown that an a^ 
monopole creates 2 units of A^ gauge charge, and vice 
versa. Thus such monopole events simply mean that 
fafl an d w a w /3 pairing terms exist in the effective ac- 
tion. These pairing terms, or monopole events, break 
the U(l) gauge symmetry down to Z 2 . 



then the most general 0(4) transformation leaving 
= 1 is U -> VlUVr, where Vl and Vr 



wi 



\w 2 \ 



are both SU{2) rotations (Vr, is spin rotation), and 
0(4) ~ SU(2) L x SU(2) R . In this representation, the 
vectors ni, ri2, h are the 1st, 2nd and 3rd columns of a 3 
by 3 rotation matrix R: 18 



R ab = TT ^ a a U(J b^ 



We should focus on one set of dual variables /(V^) 
and w(VJ) to write down the effective theory. First take 
a look at the V A operator: what is the symmetry of 
the corresponding /-spinon pairing term? The quantum 
number of this monopole operator is determined by the 
WZW term: a Skyrmion of n carries an s-wave pair of 
/-spinon. Therefore we showed that h disordered phase 
is nothing but the SPS, whose mean-field ansatz is given 
in Eq.©. 

The phase transition from SPS to CAF phase is de- 
scribed by w-boson condensation in Eq.([S]). But due 
to the topological Berry's phase, this is a rather "high- 
energy" description of the criticality. To find a low energy 
description without topological terms, we should resort 
to another formulation as shown in the next section. 

Duality between Schwinger-fermion and 
Schwinger-boson representations 

In this section we focus on the dual variables of /- 
spinons: the w-bosons. The SPS phase is then a Zi phase 
with w-bosons as Zi charges, but /-spinons as visons. In 
this formulation SPS-CAF phase transition is naturally 
presented as a Higgs condensation of w-bosons. 

First we need to represent the order parameters of the 
CAF phase in terms of w. The QSH order is h = w^aw, Under duality transformation 



(10) 



Clearly, to leave R invariant, the transformations Vl,r 
must be ±1. 

These symmetry transformation rules allow us to re- 
veal the connection between the SPS state here and the 
0-flux state in the Schwinger-boson representation ob- 
tained by Wang 3 . In Wang's work, the Neel order is 
represented by the z-boson as N = z^az in the effective 
theory. From Eq . (TTU|) . we can easily construct the du- 
ality transformation between the w-boson and z-boson 
representations: U w = U z Vr, Vr = e l ^ (7y , namely: 

w = —i={ z — i<J v z*) or z = —=(w + ia y w*) (11) 
V 2 y2 



h\ + ih 2 = (icr y w*)l e <7 ,pWp. 



but what is the Neel order parameter? Neel order in 
CAF phase corresponds to the monopole of a^, namely a 
pairing of w-boson. There are two spin-1 bosonic pairing 
order parameters satisfying this requirement, i.e. the real 
and imaginary part of (ia y w*Y aw: 

(7) 

It is easy to verify that h = hi x n 2 , so there are only 
two independent vectorial order parameters. The issue 
is, which one is the Neel order parameter N: h\ or n.2? 

A 17(1) gauge transformation w — > e l6 w generates a 
rotation in the hi,U2 plane. By fixing a proper gauge, 
we can always choose hi as the Neel order. We will work 
within this gauge N = hi throughout the phase transi- 
tion. Such a gauge fixing breaks the 17(1) gauge redun- 
dancy down to Z 2 : w — > ±w. 

The physical symmetries of the QSH (or vector spin 
chirality) and the Neel order parameters completely de- 
termine the transformation rules of the w-boson up to a 
Zi gauge redundancy: 



N = Ke[(ia y w*y aw] — z'az, 



h = w* aw 



-Re[(ia y z*yaz] 



(12) 



From Eq.([8|). (|lll) . we can obtain transformation rules of 
z-bosons: 



Ti,T 2 
T 
a 

Ca 



z — » z, 
z — > a y z, 
z —$■ ia y z* 
z —$■ a v z* . 



(13) 



Ti,T 2 


h — > h, 


N -> N, 


w — > w, 


T 


h — > h, 


N -> -N, 


w — > iw*, 


a 


h — > —h, 


N -> -N, 


w — > ia y w 


c e 


h — >• h. 


N -> -N 


w — > iw. 



(8) 

where time-reversal transformation T is anti-unitary. 
The reason why there are no further arbitrariness on the 
transformation rules of w can be easily understood by the 
following construction. If we write w-boson as an SU(2) 
matrix: 



U 



«'2 



"'9 



(9) 



which are exactly the transformation rules found by 
Wang 3 for 0-flux state up to a Z 2 gauge arbitrariness. 
This explicitly confirms that the z-bosons constructed in 
Eci. tfTTj) are the same z-bosons discussed by Wang, and 
the SPS phase here is identical to the 0-flux phase in 
Schwinger-boson description. 

Following the discussion in Ref. y, we can write down 
the general symmetry-allowed effective theory for the 
phase transition in terms of z-boson: 

L HcU| 2 + c 2 |Vz| 2 + m 2 |z| 2 + w(|z| 2 ) 2 

+ \ H (ia v z*Y [Y,(di ■ V) 3 ] z + h.c. (14) 

i 

For instance, the single time derivative term z^d T z is for- 
bidden by cr, and z T (—ia y )d T z is forbidden by Cq. Here 
Xh is the Higgs coupling which reduces the U{\) gauge 
degrees of freedom in the z-boson formulation down to 
Z 2 - d\ — —Si,d,2 — 3,2, and g^ = S1—S2 as shown in FIG. 
[5J The Higgs coupling can also be written as a pairing of 
w-bosons: A# (ia y w*y [^,-(dj • V) 3 ] w + h.c. By naive 



power counting Xh is irrelevant, therefore we have an 
0(4) critical point between the CAF (z-condensed) phase 
and the SPS (z-gapped) phase. The critical behavior of 
this transition is well-studiecU2r— . 

Discussion In this study, our main prediction is the 
CAF phase. Unlike the usual AF phase, CAF phase has 
two order parameters: Neel N and QSH h. As CAF 
phase is very likely to be the magnetic ordered phase 
adjacent to the spin liquid phase found in the Hubbard 
model on the honeycomb lattice, in the following we pro- 
pose an explicit numerical methods to detect the CAF 
phase. 

One can directly measure the QSH order by (n(x) ■ 
n(0)) correlation function, or the vectorial spin chirality 
correlation function {(v i+x j +x S i+x x S j+X ) ■ [v i:j Si x Sj)). 
Because QSH order is odd under a ■ T, while pure Neel 
order is a • T even, one does not expect a long range cor- 
relation of QSH order in a usual AF phase. Thus the long 



range QSH correlation function is an intrinsic signature 
of the CAF phase. In addition, one can show that the 
QSH direction is normal to the Neel direction. For ex- 
ample, one can pin the Neel order by an infinitesimal (in 
thermodynamic limit) staggered magnetic field along S z 
direction, and then to measure the QSH order parameters 
by correlation function. One should find the QSH order 
parameters have only x, y components. In the real world, 
as mentioned in Ref. |2|, such an exotic spin liquid may 
be realized in many candidate systems: e.g. expanded 
graphene-like system in group IV element o 23 ' 24 , as well 
as fermions in optical lattice o 25 ' 26 . 

YR thanks for helpful comments from Ashvin Vish- 
wanath. YML thanks Prof. Ziqiang Wang for sup- 
port during this work under DOE Grant DE-FG02- 
99ER45747. YR is supported by the startup fund at 
Boston College. 
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Appendix A: Derivation of the mutual 
Chern-Simons term 



We start from the following low energy effective La- 
grangian of spinon fields ip (see Eq. [3]) in imaginary time 
(i.e. Euclidean space-time): 



C 



eff 



ip[hn(dii. - i<v) + im n ■ a ]tJj (Al) 



where we define i/S = 7/^7° . For simplicity let's de- 
note — iQ^ 1 = 7^(<9 M — ia M ) + ran • a, then integrat- 
ing out spinon fields if) yield the effective action S = 
-lndet(^- 1 ) = -Trln(^- 1 ). Following the spirit of 
Abanov and Wiegman n 16 ' 27 , we use large-m expansion 
to obtain the low energy effective theory in the longwave- 
length limit lj <C m. By defining Qq — i (7 M <9 M + mh ■ a) 
we have Q~ l = Q$ + ^V7^- Let's denote ^ = 7^9^ and 
similarly <ji = j^a^ and we have 



S 



-Trln(g - 1 + 
= 5 



) = -Tr 111(0 -Tr\n(l + g o9 
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Here So gives the nonlinear-sigma-model dynamics ' 
(d^n) 2 of vector n, while the coupling between vector 



and gauge field a M is given by the 2nd term. In the large- 
m expansion we consider only the leading-order term: 

Si = -TrQSot) = -T*{Go 1 [Go 1 (Go l )T 1 <f} 

It's straightforward to check that Gq 1 ^^ 1 )^ = —d 2 + 
m 2 — ma ■ c/)h, therefore large-m expansion leads to 
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-d 2 + m 2 ' 



and consequently 
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i + mh ■ (?) / ma ■ q)h .. i 
-d 2 + m 2 —d 2 + m 2 



It turns out that I = 0, 1 terms both vanish and the 
leading-order correction to the low energy effective action 
is the following topological term: 



Notice that in the CP 1 parametrization of order pa- 
rameter n — w^aw, spinor w is the eigenvector of h ■ a 
whose spin orientation is along unit vector n. Therefore 
the Skyrmion current of n 
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w\n ■ (d v n x d\h) 



(A3) 



which equals half the winding number of n wrapping 
around S 2 , is nothing but the Berry's phase 28 for spinor 
w. Since spinor w obtains 7r phase (i.e. a minus sign) as 
h wraps around S 2 once (i.e. h covers 4ir solid angle), 
this gives a direct correspondence between the Skyrmion 
current density and the U (1) gauge field strength F^ v 
coupled to w: 
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\h ■ (d v n x d\h) 



^^u\d v Ax 



(A4) 



Therefore the topological term (|A2|) is exactly the mutual 
Chern-Simons term mentioned in ([3]) 
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